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Abstract 

We report results of density functional theory (DFT) calculations on some colorless aromatic systems adsorbed on a Ti02 
nanocluster, in order to explain experimental results regarding the photocatalytic degradation of these pollutants under visible light 
irradiation. Based on our modeling, we are able to clarify why transparent pollutants can degrade under visible light in the presence 
of a catalyst that absorbs only in the UV, to explain experimental data regarding differences in the efficiency of the degradation 
process, and to state the key requirements for effective water-cleaning. For that purpose, we analyze the absorption spectrum of the 
free and adsorbed molecules, the binding configurations, the matching of the energy levels with the oxide catalyst and the likeli- 
hood of the charge-transfer to the substrate. The comparison between several colorless aniline and phenolic systems allows a corre- 
lation between the chemical structure and the degradation rate of these pollutants. 



Introduction 

Titania, TiC>2, has been widely used as photocatalyst for envi- 
ronmental applications [1-6], particularly for waste water purifi- 
cation. Due to its large band gap Ti02 absorbs only UV radia- 
tion, a fact that limits the efficiency and keeps the costs of the 
photocatalytic degradation of environmental pollutants high. To 



be used under visible light irradiation, in the range of wave- 
lengths where the solar spectrum has its maximum, the elec- 
tronic band structure of the photocatalyst has to be modified in 
various ways [6]. Alternative approaches to the modification of 
the TiC>2 photocatalyst are the self-sensitized degradation of 
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dyes which absorb visible light [7,8] and the photocatalytic deg- 
radation of colorless organic compounds by formation of a 
charge-transfer-complex, CTC [9,10]. The assumption of a 
surface CTC in the visible light catalysis was supported by 
subsequent work on various other types of systems, such as 
phenolic compounds [11,12], fluoroquinolone antibacterial 
agents [13], and various colorless aromatic pollutants [14]. 

Despite the extensive experimental work, the role of the key 
factors that influence the microscopic mechanism of the photo- 
catalytic processes is not entirely understood. Our goal is to 
answer a few key questions, regarding photocatalytic degrada- 
tion: i) Why can transparent pollutants degrade under visible 
light in the presence of a catalyst that absorbs only in the UV? 
ii) Why are some pollutants degraded more efficiently than 
others? iii) What are the requirements for an effective water- 
cleaning process? To answer these questions we start from an 
analogy with the photoelectrochemical Gratzel cells [15,16]. 
We argue here that efficient photocatalytic degradation of pollu- 
tants under visible light irradiation has to meet similar require- 
ments to the ones of the dyes in Gratzel cells. In particular, the 
anchoring mode of the pollutant to the Ti02 surface influences 
the electron transfer [17]. The most commonly used anchoring 
group is the carboxylic acid group (-COOH) [16]. It ensures 
strong binding of the dye on the surface and promotes the 
charge transfer. The anchoring of the salicylate group on Ti02 
has also been studied [18,19]. At the surface, both substituent 
groups of a benzene derivative are involved in the complexa- 
tion of colloidal titanium dioxide [18]. This results in the forma- 
tion of a six-atom ring with a chelating type of bonding to the 
same Ti(IV) ion. Similarly, the binding of the salicylic acid to 
titania was thought as bidentate chelate through the oxygen 
atoms of -OH and of-OCOH [14,20]. 

Theoretically, density functional theory (DFT) calculations 
showed [21-23] that the binding of the carboxy group to titania 
is bidentate bridging, with the monodentate anchoring being 
less stable [24-27]. The higher performance of the dyes with 
both carboxy and hydroxy anchoring groups [28] has led us to 
revisit earlier studies of a dye with three types of anchoring 
groups: -OH, -COOH and -S0 3 H [29,30]. We showed that 
although the salicylate does use both the carboxy and hydroxy 
substituent groups, the binding configuration is not bidentate 
chelate, as previously thought [14,18-20]. 

Building upon the experience gained while modeling materials 
for photoelectrochemical cells, we report here results of DFT 
and time dependent DFT (TD-DFT) calculations performed on 
several colorless aromatic pollutants, as well as complex 
systems consisting of benzene derivatives adsorbed on a Ti02 
nanocluster. To answer the questions raised above we deter- 



mine the electronic structure and the optical spectra of the 
pollutant itself, and find where the deprotonation is more likely 
to take place. We also simulate the pollutant-catalyst system to 
analyze the binding configurations. We discuss the energy level 
alignment between the pollutant and the catalyst as well as the 
charge transfer between the pollutants and the oxide. We 
compare our theoretical results with the experimental data avail- 
able, particularly with the work of Wang et al. [14] on phenol 
(Ph), benzoic acid (BA), /?-hydroxybenzoic acid (pHBA) and 
salicylic acid (SA), which attempted a correlation between the 
efficiency of photocatalytic degradation and the chemical struc- 
ture of the pollutants over Ti02- 

Results and Discussion 

This section is divided in five parts. The first describes the 
computational details whereas the second focuses on the opti- 
mized geometry and electronic structure of the free pollutants. 
The third subsection presents the binding of the pollutants to the 
titania nanocluster, the fourth presents the optical properties of 
the adsorbed pollutants, and the last subsection attempts to 
explain the experimental data as well as the key requirements 
for efficient photodegradation based on theoretical arguments. 

Computational details 

The structures of all pollutants were optimized in neutral as well 
as deprotonated forms, using DFT [31-33], with the B3LYP 
exchange-correlation functional [34,35] and the double-^ 
DZVP basis set including polarization functions for the valence 
electrons [36] of the free pollutants. In the case of the more 
complex pollutant-catalyst system the less demanding 3-21G(d) 
basis set was used for geometry optimization. All optimized 
structures were checked for stability by means of vibrational 
analyses. Time-dependent DFT [37] calculations of the molecu- 
lar orbitals and the electronic transitions were performed in 
water by means of the polarizable continuum model (PCM) 
[38,39]. We used the same B3LYP functional and TZVP basis 
sets [36]. In the case of the pollutants adsorbed on the catalyst, 
the electronic states were accurately computed by using DZVP 
basis sets [40]. The Gaussian03 package [41] was used in all 
calculations. 

Free pollutants - electronic structure and 
optical properties 

During the photocatalytic degradation the benzene derivatives 
(phenol, Ph, benzoic acid, BA, /?-hydroxybenzoic acid, pHBA, 
and salicylic acid, SA) undergo various processes, including 
deprotonation. We first perform a comparative analysis of the 
likelihood of the deprotonation process. For that, we found the 
optimized geometries and determined the total energy for all 
four aromatic pollutants in neutral and their deprotonated forms. 
The simulations of the deprotonated forms were performed by 
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taking away a proton from the anchoring group, as it generally 
happens when the molecule is bound to the substrate. The opti- 
mized structures are shown in Figure 1 . 

Obvious differences between the geometries of the neutral and 
deprotonated forms are observed for salicylic acid. In the case 
of doubly deprotonated SA the carboxy group moves away 
from the plane of the aromatic ring by 86.2°, whereas for the 
form that is singly deprotonated at the carboxy group the torsion 
angle is 48.7°. A more careful analysis also reveals some slight 
differences between the bond lengths and some small distor- 
tions of the bond angles but these aspects are not of crucial 
importance in the following. The key question is at which 



anchoring group the deprotonation is more likely to take place. 
To find the answer we can look at the total energy of the pollu- 
tants as well as at the proton affinity, PA, equal to the differ- 
ence between the energy of the deprotonated form and the 
energy of the neutral compound [42], reported in Table 1. For 
both pollutants, the lowest energy between the two singly 
deprotonated forms is obtained when the hydrogen atom of the 
carboxy group is removed. 

The DFT calculated energies of the key molecular orbitals of 
the pollutants in neutral and deprotonated forms are repre- 
sented in the diagram shown in Figure 2. To put everything in 
perspective we also represented the results of DFT calculations 




Figure 1: Optimized structure of pollutants in neutral and deprotonated forms calculated at DFT/B3LYP/DZVP level: a) neutral and b) deprotonated 
Ph, c) neutral and d) deprotonated BA, e) neutral, f) and g) deprotonated, h) doubly deprotonated pHBA, i) neutral, j) and k) deprotonated, I) doubly 
deprotonated SA. 
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Table 1 : Proton affinities of pHBA and SA in neutral and various 
deprotonated forms, in water, based on DFT calculations at B3LYP/ 
TZVP level. The labeling corresponds to the one used in Figure 1. 



pHBA 


(e) 


(f) 


(g) 


(h) 


PA (eV) 




12.97 


12.81 


26.00 


SA 


(i) 


G) 


(k) 


(I) 


PA (eV) 




13.01 


12.78 


26.09 



for Ti02 clusters [23], particularly the edges of the valence and 
conduction bands of titania. As discussed in [23,43], and the 
references therein, DFT calculations on finite size clusters tend 
to overestimate the gap with respect to the experimental values. 
It can be seen that the energy difference between the lowest 
unoccupied molecular orbital (LUMO) and the highest occu- 
pied molecular orbital (HOMO) is larger for all pollutants than 
the calculated band gap of the semiconductor. As a result, the 
absorption spectra of the free pollutants have peaks further in 
the UV region than the catalyst. Another important observation 
is related to the energy level alignment between the pollutants 
and the catalyst. As all LUMOs of the pollutants are above the 
conduction band edge of Ti02, the charge transfer to the semi- 
conductor is possible. We note that the energies of the deproto- 
nated forms are shifted upwards with respect to the neutral 
pollutants, as expected due to the missing proton. 

The electronic spectra of the pollutants, simulated by TD-DFT 
calculations, are displayed in Figure 3. The spectra are all in the 
UV, as expected since the pollutants are all colorless, and in 
agreement with the results of DFT calculations presented in 
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Figure 3: Simulated UV-vis absorption spectra of the first 10 transi- 
tions for neutral pollutants in water, calculated by TD-DFT. The spec- 
tral lines were convoluted with Gaussian distributions of 20 nm 
linewidth at half maximum. 



Figure 2. The energy of the electronic transitions, obtained by 
TD-DFT is a better indicator of the gap than the 
LUMO-HOMO energy difference determined by DFT [23]. 

Of all pollutants, salicylic acid has peaks at higher wavelengths, 
just below 300 nm. The spectrum of Ph has a HOMO^LUMO 
transition at 248 nm, a HOMO->LUMO+l transition at 215 nm 
and HOMO-1— >LUMO transitions at 187 nm. In the case of 
BA the transitions are at 259 nm (HOMO— >LUMO), 237 nm 
(HOMO-l^LUMO), and 192 nm (HOMO^LUMO+1). 
Similarly, for pHBA we found transitions at 258 nm 
(HOMO^LUMO), 217 nm (HOMO^LUMO+1), and 191 nm 
(HOMO-l^LUMO+1). Finally, in the case of SA the spec- 
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Figure 2: Diagram of energy with respect to vacuum as well as potential versus the normal hydrogen electrode (NHE) of pollutants in both neutral and 
deprotonated forms, calculated at DFT/B3LYP/TZVP level in solution. In the case of pHBA and SA the deprotonation is on the carboxy group. 
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trum is closest to the visible range: 289 nm (HOMO— >LUMO), 
234 nm (HOMO-ULUMO), and 207 nm 
(HOMO— >LUMO+l). All transitions have 71-71* character, for 
all pollutants. The HOMO— >LUMO transitions are weak, the 
most intense peaks being located for all four compounds in the 
range of 187-207 nm. These strong transitions involve levels 
just below HOMO and/or just above LUMO. 

Before ending this section we present in Figure 4 the key mo- 
lecular orbitals of the pollutants, as they will be useful in identi- 
fying in a later section the MOs of the complex pollutant-cata- 
lyst system. 

Pollutant adsorption - binding configurations 

As argued in the Introduction, the anchoring modes of the 
pollutants to the catalyst surface are of crucial importance for 
the charge transfer. But before discussing adsorption, a few 
comments regarding the titania cluster are in place. We model 
Ti02 nanoparticles by a geometry optimized cluster with the 
molecular formula ^2405(^4. Prior to optimization, the cluster 
was cut from the experimental anatase structure with (101) and 



(001) surfaces [44]. The optimization led to some slight distor- 
tions from the lattice geometry, which lowered the surface 
energy. To prevent the occurrence of dangling bonds for some 
peripheric oxygen atoms and to avoid the problem of the 
surface states in the gap [45] we introduced four hydrogen 
atoms [29,46]. The resulting Ti2405oH4 cluster has a length of 
12.76 A and a width of 7.39 A and provides a reasonable 
compromise between accuracy and computational costs given 
the small size of the pollutant molecules [29]. One last 
comment on this topic is related to the nature of the states in the 
valence and conduction band of titania. The valence band is 
dominated by the contributions from the p orbitals of oxygen, 
whereas the conduction band is dominated by the d orbitals of 
the Ti atoms [29,46]. 

We start our discussion of adsorption with Ph, which, after 
deprotonation, can bind to the (101) surface of the substrate 
through the oxygen atom. This is shown in Figure 5, which 
displays the optimized structure resulting from the DFT calcula- 
tions. We note that the geometry optimization was performed 
for all pollutants starting with the molecule distanced from the 
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Figure 4: Isodensity surfaces (0.03 e/bohr 3 ) of the key molecular orbitals of the deprotonated pollutants, calculated at DFT/B3LYP/ TZVP level in 
water. In the case of pHBA and SA the deprotonation is on the carboxy group. 
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pHBA (h) 




pHBA (c) 



Figure 5: Optimized geometry of the pollutant-catalyst systems calculated at DFT/B3LYP/3-21G* level; pHBA is anchored either through the hydroxy 
(h) or the carboxy (c) groups. 



cluster. The final geometry with pollutants bound to the Ti 
atoms of the (101) surface was robust with respect to variations 
in the initial configuration. The distances and angles relevant to 
the binding of Ph to the titania cluster are reported in Table 2. 
The dihedral angle x, which is measured between the plane of 
the ring and the surface of the cluster (defined by three adjacent 
Ti ions), is 71.1° tilted away from the normal plane. 

In the case of BA, the adsorption is through the carboxy group. 
Experimental studies have shown that carboxylic acid groups 
can have various binding configurations with the Ti(IV) ions, 



ranging from monodentate ester-like binding [47], bidentate 
bridging [48], or both bidentate chelate and bridging [49,50]. 
However, theoretical calculations demonstrated that the 
preferred anchoring is bidentate bridging, with one proton trans- 
ferred to a nearby surface oxygen [22,24-26]. Our calculations 
for BA are consistent with these earlier results, the bidentate 
bridging being the preferred adsorption mode, as shown in 
Figure 5b. The Ti-0 bond distances reported in Table 2, are 
different, 2.058 A and 1.991 A the axis of the molecule being 
slightly tilted. The dihedral angle is 86.1° almost perpendicular 
to the surface of the catalyst. 



Table 2: Bond distances and bond angles relevant to the binding of the pollutants to the catalyst, after geometry optimization at DFT/B3LYP/3-21G* 
level; pHBA is anchored either through the hydroxy (h) or the carboxy (c) groups. 



parameter 



Ph 



BA 



pHBA (h) 



pHBA (c) 



r(Ti-0) 
r(Ti-01) 
r(Ti-02) 
e(Ti-O-C) 
9(Ti-01-C) 
6(Ti-02-C) 



1.835 A 



145.5° 



71.1° 



2.058 A 
1.991 A 

128.3° 
131.0° 
86.1° 



1 .940 A 



120.2° 



13.6° 



2.049 A 
1.985 A 

128.7° 
131.0° 
84.0° 
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For pHBA there are two possible anchors, one is the carboxy 
group, as in the case of BA, the other is the hydroxy group as 
for Ph. Our simulations provided different total energies for the 
two cases, the preferred binding being through the carboxy 
group, by an energy difference of about 1 eV. Such a result is 
not a surprise, as the strength of the mechanical adhesion is 
higher when two Ti-0 bonds are involved instead of just one 
[45]. We also note the different orientation of the aromatic ring 
in the two cases. When the binding is through the carboxy 
group the aromatic ring is almost perpendicular to the surface 
of the catalyst (x = 84.0°). In contrast, when the binding is of 
the (h)-type the ring is tilted almost parallel to the surface 
(x = 13.6°), reminding of n stacking interactions. 

In the case of SA adsorbed on T1O2 previous studies claimed 
that in the salicylate complex formed at the interface both 
substituent groups are involved, which leads to the formation of 
a six-atom ring with a chelating type of bonding to the same 
titanium ion [18,19]. Similarly, the binding of SA was thought 
as bidentate chelate through the oxygen atoms of -OH and of 
-COOH [14,20]. In a previous paper [29] we showed that for 
Mordant Yellow- 1 0 anchoring through the S A part can lead to 
three binding modes, depending on the degree of deprotonation, 
but none of them was bidentate chelate. Returning to the present 
study, for SA we performed geometry optimization calcula- 



tions in three cases. We deprotonated the carboxy group, case 
(c), the hydroxy group, case (h), and both, case (c&h), when we 
left one proton to find its equilibrium position on the Ti02 
surface. All these cases are illustrated in Figure 6 and Table 3. 

The first question to answer is which of the three types of 
anchoring is most stable. Our DFT calculations showed the 
most stable configuration is (c&h). At about 0.58 eV higher lies 
the SA deprotonated at the hydroxy group (h), followed at 
1.18 eV by the configuration obtained by deprotonation of the 
carboxy group (c). Consequently, we confirm our previous 
calculation on Mordant Yellow 1 0, showing that the deprotona- 
tion of both groups lowers the energy the most, allowing for 
stable 'mechanical' structures with three pillars of unequal bond 
lengths to three different but adjacent Ti(IV) ions. The shortest 
bond length is the one to the oxygen of the hydroxy group. The 
triple binding is allowed by a rotation of x' = 62.0° of the 
-COO~ group with respect to the aromatic plane. On the other 
hand, the benzene ring encloses a dihedral angle of 72.0° with 
the plane of the surface. 

Next, when the deprotonation takes place at the hydroxy group 
alone, a second bond is formed through the oxygen atom of the 
carboxy group and the third is a hydrogen bond involving the 
same carboxy group. We clearly demonstrate that the 'common 




Figure 6: Optimized geometry of SA adsorbed on titania, calculated at DFT/B3LYP/3-21G* level. 



(c&h) 



Table 3: Distances and angles relevant to the binding of SA to the catalyst, after geometry optimization at DFT/B3LYP/3-21G* level. 



parameter 



SA (c) 



SA (h) 



SA (c&h) 



r(T\-0) 
r(Ti-01) 
r(Ti-02) 
r(0-H) 
e(Ti-O-C) 
9(Ti-01-C) 
9(Ti-02-C) 

T 
T' 



2.051 A 
1.969 A 



128.1° 
133.0° 
84.8° 



1.850 A 

2.031 A 
1.112 A/1.386 A 
146.0° 
123.6° 

45.0° 
24.6° 



1.830 A 
2.017 A 
2.036 A 

163.2° 
121.2° 
125.0° 
72.0° 
62.0° 
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knowledge' that the binding configuration is in a six ring 
chelate, suggested by infrared spectroscopy measurements 
[19,20] or simply assumed by other authors [14,18], is wrong. 
There are two direct bonds to two different Ti(IV) ions and one 
indirect bond to an oxygen through the H atom. Again, the third 
bond is allowed by a rotation of the -COO~ group by t' = 24.6°. 
Overall, the plane of the ring is tilted by 45.0° with respect to 
the surface of the catalyst. The last case involves only the usual 
bidentate bridging through the carboxy group. This anchoring 
configuration is overall similar to what we observed for BA and 
pHBA (c), with comparable bond lengths, bond angles and 
dihedral angles. 

Adsorbed pollutants - optical properties 

We recall our conclusion drawn based on Figure 2 and Figure 3, 
that the absorption bands of the free pollutants are situated 
further into the UV region than those of the catalyst. In contrast, 
the TD-DFT simulated optical spectra of the pollutants 
adsorbed onto the catalyst, displayed in Figure 7, show for all 
compounds a strong shift toward higher wavelengths. This red- 
shift is what makes possible the photocatalysis under visible 
light irradiation. In order to better understand why this red-shift 
takes place, it is useful to plot the densities of states for all 
adsorbed pollutants. But before that we take a more careful look 
at Figure 7. 
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Figure 7: Simulated UV-vis absorption spectra of the pollutants bound 


to the TiC>2 cluster, calculated by TD-DFT in water. The spectral lines 


were convoluted with Gaussian distributions of 20 nm linewidth at half 


maximum. 







The first (low energy/high wavelength) bands correspond for all 
systems to HOMO— ►LUMO transitions. However, these transi- 
tions have very small oscillator strengths, as it can be seen in 
Table 4. For instance, the only adsorbed pollutant with a band 
beyond the blue region of the spectrum is Ph, whose first tran- 
sition is at 524 nm, but the intensity is extremely weak. The 
stronger UV-vis absorption bands of the adsorbed pollutants are 



Table 4: Wavelength, oscillator strength and character of the most intense optical transition for the four pollutants adsorbed on the catalyst, calcu- 
lated at TD-DFT/B3LYP/DZVP level in water. 

pollutant A (nm) f character 



524 0.0002 HOMO^LUMO 

445 0.034 HOMO^LUMO+2 

Ph/T l24 0 5 oH4 387 0Q38 HOMO-LUMO+13 

345 0.029 HOMO^LUMO+36 



354 0.001 HOMO^LUMO 

348 0.019 HOMO-1^LUMO 

BA/Ti 24 0 5 oH4 33Q Q115 HOMO-1-LUMO+2 

306 0.041 HOMO-1^LUMO+4 



pHBA/Ti 24 0 5 oH 4 (c) 



425 
396 
365 
335 
303 



0.013 
0.114 
0.040 
0.055 
0.018 



HOMO^LUMO 
HOMO^LUMO+2 
HOMO^LUMO+6 
HOMO^LUMO+13 
HOMO^LUMO+29 



SA/Ti 24 O 50 H 4 (c&h) 



454 
428 
399 
346 
319 



0.002 
0.004 
0.014 
0.047 
0.008 



HOMO^LUMO 
HOMO^LUMO+2 
HOMO^LUMO+8 
HOMO^LUMO+19 
HOMO-1^LUMO+6 
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located below 475 nm. Except for BA, whose key orbital seems 
to be HOMO-1, for the other adsorbed pollutants the transitions 
with noticeable oscillator strengths are from the HOMO to 
states above the conduction band edge, such as LUMO+2. 

For the adsorbed Ph we observe three major peaks, one in the 
visible range, at 445 nm, the other two in the UV range. In 
the case of the pHBA on the Ti02 nanocluster, the 
HOMO— >LUMO band at 425 nm has a significant intensity, 
still much smaller than the next. At the limit of the visible range 
there is a transition at 396 nm, from HOMO to LUMO+2. 
When bound to the titania cluster, SA has two weak bands in 
the visible range and stronger absorption in the UV. Finally, the 
adsorbed BA has only absorption bands in the UV. 



In order to better understand the electronic spectra we looked at 
the electronic density of the key molecular orbitals, displayed in 
Figure 8. As expected, the HOMOs have most of the charge 
located on the pollutant, whereas the LUMOs correspond to the 
conduction band edge of the catalyst. The excited state of the 
adsorbed pollutant, identified by comparison with the electron 
densities of the free pollutants represented in Figure 4, is situ- 
ated deep into the conduction band for all compounds. For 
instance, in the case of BA/Ti2405oH4 the corresponding MO is 
the 52nd from the LUMO, whereas in the other cases there are 
more than 80 states between the conduction band edge of titania 
and the excited state of the pollutant. Of crucial importance in 
the absorption spectra are the MOs mentioned in the transitions 
reported in Table 4. A common feature of these orbitals is the 
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Figure 8: Isodensity surfaces (0.03 e/bohr 3 ) of the key molecular orbitals of the deprotonated and adsorbed forms of pollutants, calculated at DFTV 
B3LYP/DZVP level in water. 
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mixed character, the electron density extending from the semi- 
conductor onto the pollutant, as it can be seen from the second 
row of Figure 8. In fact, it is this derealization of the charge 
from the pollutant onto the catalyst in the case of the HOMOs 
as well as from the semiconductor to the adsorbed molecule, in 
the case of higher empty states, that make possible the higher 
wavelength electronic transitions and the absorption in the 
visible range. Although taken separately both the pollutants and 
the catalyst absorb in the UV, the states with mixed character 
lead to allowed optical transitions. 

We can even better understand these concepts looking at the 
density of states (DOS) for the adsorbed pollutants presented in 
Figure 9. For all four systems, the densities of states have some 



common features. For instance, the edge of the p-type valence 
band is located just above -7.9 eV, whereas the d-type conduc- 
tion band edge is situated just below -3.7 eV. We note that the 
DFT calculated gap of about 4.2 eV for TiC>2 is overestimated 
with respect to the experimental value [23,43]. 

Common to all compounds is that in the gap, significantly 
higher than the valence band edge, there are two occupied 
states, HOMO-1 and HOMO, well localized on the pollutants. 
These two orbitals are separated by about 0.7 eV for three of the 
four systems, the only exception being BA, for which the two 
states are almost degenerate (the energy difference is only 
0.08 eV). For the three systems for which the HOMO-1 is 
significantly below the HOMO, only the HOMO plays an 
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Figure 9: Density of states of the pollutant-catalyst complex systems, calculated at the DFT/B3LYP/DZVP level in water. Bottom scale is for the 
contribution of the Ti2405oH 4 nanocluster (dotted line) whereas the top scale for the contribution of the pollutant (continues line). For all systems, the 
edge of the p-type valence band of the semiconductor is located at about -7.9 eV, whereas the d-type conduction band edge is situated at about 
-3.7 eV. The HOMOs of the pollutant are located in the gap of the semiconductor. 
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Table 5: Contributions of pollutant, anchor group, and substrate, in %, 



to the electron density of MOs involved in lowest electronic transitions, 


caicuiatea at ui- 


/boLYr/UZ.Vr level. 






MO 


pollutant 


anchor 


Ti 2 405oH 4 


phenol 


HOMO 


86.95 


19.09 


13.05 


LUMO 


0.14 


0.02 


99.86 


LUMO+2 


1.29 


0.63 


98.71 


LUMO+89 


42.65 


0.01 


57.35 


benzoic acid 


HOMO 


99.58 


0.67 


0.42 


LUMO 


0.57 


0.38 


99.43 


LUMO+2 


0.15 


0.01 


99.85 


LUMO+52 


11.41 


4.61 


88.59 


p-hydroxybenzoic acid 


HOMO 


96.29 


8.86 


3.71 


LUMO 


0.45 


0.29 


99.55 


LUMO+2 


0.04 


0.001 


99.96 


LUMO+81 


91.59 


0.32 


8.41 


salicylic acid 


HOMO 


88.33 


21.84 


11.67 


LUMO 


0.27 


0.14 


99.73 


LUMO+8 


2.48 


1.49 


97.52 


LUMO+84 


31.81 


2.38 


68.18 



important role in the optical spectra. In the fourth case, both 
HOMO-1 and HOMO can play a role but lying lower in 
energy, closer to the valence band edge, the energy difference is 
high and the transitions are in the UV. The energies of the 
HOMOs are: -6.36 eV, -7.39 eV, -6.78 eV and -6.68 eV for 
the complex systems with adsorbed Ph, BA, pHBA, and SA, 
respectively. 

In line with the previous discussion, important to note is the 
mixed character of the states involved in the optical transitions. 
Looking carefully at the density of states, we can observe that 
the HOMOs have some small contribution from the catalyst 
and, similarly, the key states in the conduction band have a 
small contribution from the pollutant. The weak DOS peaks in 
the conduction band (note the different scale used) correspond 
to states with 7t* character and with sizeable electron density 
distributed over the pollutant, as shown in Figure 8. 

For a more quantitative image we report in Table 5 the contri- 
butions to the electron density of the MOs involved in the 
optical transitions. This way we can see that the charge of the 
HOMOs is not entirely distributed on the pollutant. While BA 
keeps over 99% of the charge, the other pass more than 13%, 
3%, and 11% (for Ph, pHBA, and SA, respectively) to titania. In 
the conduction band, the picture reverses, as most of the charge 
is on titania. However, even in this case there are states with a 
sizeable charge located on the pollutant. 

The data reported in Table 5 allow us to draw additional conclu- 
sions regarding the likelihood of the electron transfer. In 
Marcus' theory of the electron transfer [51-53] an important 
factor in the expression of the electron transfer rate is the elec- 
tronic matrix element describing the electronic coupling 
between the excited state of the pollutant and a state in the 
conduction band of the catalyst. When the orbitals of the two 
separate components are known, an indicator of the strength of 
the matrix element may be the overlap between those states. 
When the orbitals of the entire system are available, as it is the 
case in our work, we can alternatively look at the degree of 
mixing indicated by the electron density present on each 
component and at the electron density on the anchoring group 
that binds the pollutant to the surface of the catalyst. Taking 
another look at Figure 8, we note that the flow of charge from 
the pollutant to the catalyst has to take place through the anchor. 
If the electron density on the anchor is high in the donor state, 
the tendency for charge transfer has to be stronger. The states 
with little density on the atoms of the anchoring group (such as 
the HOMO of pHBA/Ti02, which has a nodal plane through the 
center of the carboxy group), are less likely to favor a charge 
transfer. Table 5 reports electron densities on the anchoring 
group for each pollutant-catalyst system. Comparing the MOs 



involved in the main optical transitions, LUMO+2 for Ph and 
LUMO+8 for SA, we see that the percentage of the charge 
localized on the pollutant and on the anchoring group is about 
two times larger for SA, suggesting a more efficient charge 
transfer. 

One last comment in this section refers to the role of the 
vibronic coupling effects on the charge transfer to the Ti02 
cluster. In the case of Gratzel cells it has been shown [54,55] 
that the kinetics of the electron transfer from an excited dye to 
the titania nanoparticle may be influenced by the vibrational 
motion of nuclei. The vibronic perturbation, due to the inter- 
play of electron-electron interactions and the internal vibra- 
tions of the benzene derivative, may facilitate the charge 
transfer also in the case of photocatalytic degradation of the 
pollutants studied here. 

Discussion and comparison with experi- 
mental data 

In this section we compare our theoretical results with the 
experimental data available, particularly with the work of Wang 
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et al. [14], who showed that BA can hardly be degraded and, for 
the other three, the order of the degradation efficiency is: 
SA > Ph > pHBA. The questions we attempt to answer are: 
Why do some pollutants degrade faster than others? What are 
the requirements for an efficient photocatalytic degradation of 
the pollutants under visible light irradiation? We base our 
answer on an analogy with the photoelectrochemical Gratzel 
cells [15,16]. An efficient photocatalytic degradation under 
exposure to visible light requires: i) strong adsorption of the 
pollutant to the catalyst, ii) intense absorption of the pollutant or 
the combined pollutant-catalyst system in the visible range of 
the spectrum, iii) proper energy level alignment of the excited 
state of the pollutant and the conduction band edge of the cata- 
lyst, and iv) fast charge transfer from the pollutant to the cata- 
lyst. Other requirements would regard the chemical reactions 
that take place after the charge transfer but we will not address 
those issues here. 

Starting with the first criterion, all four pollutants studied here 
can bind to the catalyst. However, the strength of the bond is 
not the same, as at one limit SA can exhibit a triple bond, 
whereas Ph can form only a single bond. Binding through the 
carboxy group is most likely bidentate bridging, in agreement 
with various theoretical observations [22,24-26,29,46] and in 
contrast to some other opinions [56]. Also, the binding of SA 
involves both substituent groups, but does not lead to the forma- 
tion of a six-atom ring with a chelating type of bonding to the 
same Ti(IV) ion, as previously considered [14,18-20,57]. 
Instead, the proton migrates on the surface leaving all three 
oxygen atoms available for bonds to three adjacent Ti ions. 
Moreover, even if the proton were kept by the carboxy group, 
the binding would still not be chelate and an additional H-bond 
would be formed. 

Checking the four colorless pollutants studied here against the 
second criterion, we note that BA/Ti02 has extremely poor 
absorption in the visible range. The other systems have some 
absorption bands at low wavelengths, due to the mixed char- 
acter of the key MOs. For instance the HOMOs, localized 
mostly on the pollutant, have some small contribution from the 
catalyst, whereas for some states in the conduction band the 
situation is reversed. As a consequence, low energy transitions 
that were forbidden become allowed. Based strictly on the 
absorption (see Figure 7) we would be led to the conclusion that 
adsorbed Ph should degrade faster than SA (which has weak 
bands at 454 nm and 428 nm) and pHBA (425 nm). 

The third criterion, asking that the excited state of the pollutant 
lies higher than the conduction band edge of the catalyst, is met 
by all four pollutants. Here, an observation regarding a running 
controversy may be useful, especially as it paves the way for the 



discussion of the next criterion. Using arguments based on 
Marcus' theory of electron transfer and making some simpli- 
fying assumptions, it was suggested [58] that the larger the 
'driving force' (i.e., the difference in energy between the 
excited state of the molecule and the conduction band edge of 
the semiconductor) the higher the injection rate. However, there 
are some experimental observations [29,59-62] backed by some 
arguments derived through molecular modeling [29,46], which 
contradict this claim, suggesting that in practice a large driving 
force is not necessarily a guarantee for high electron transfer 
rates. 

The forth criterion, requiring a fast charge transfer from the 
pollutant to the substrate, is more difficult to analyze. One of 
the factors that influence the transfer rate is the matrix element, 
which, in turn, can be correlated with the overlap integral 
between the two states. The emphasis on the orbital overlap was 
underlined long ago [63] suggesting that the %* orbitals of the 
carboxy group would promote rapid electron injection into the 
conduction band of Ti02, due to its d-symmetry, but not that of 
Sn02 or ZnO, which have predominantly s-character. Keeping 
in mind that for Gratzel cells the charge transfer is optimized in 
the case of a strong overlap between the dye and the semicon- 
ductor [18], which is favored by the presence of charge on the 
anchor, by analogy, in the case of the pollutant-catalyst system, 
the electron density on the binding group can provide some 
information regarding the likelihood of the charge flow. We 
found that such a comparison favors SA against Ph and both 
against pHBA. 

A summary of the present discussion is displayed schemati- 
cally in Figure 10. The left panel shows the density of states and 
reveals the mixed character of the key orbitals, such that, due to 
charge derealization optical transitions in the visible become 
possible. The diagram in the central panel illustrates the energy 
level alignment and the position of the MOs involved in the 
excitation by absorption of visible light and in the electron 
transfer from the pollutant to the catalyst. Although typically 
the 'driving force' for the charge transfer is defined [53,58] as 
the difference between the energies of the excited state 
(LUMO+84, in our case) and of the conduction band edge 
(LUMO), here it may be more meaningful to consider the 
LUMO+8 orbital, which is much closer in energy and has a 
higher overlap with the LUMO. In the right panel the corres- 
ponding molecular orbitals are grouped to point out the electro- 
optical processes but they also presented the unusual anchoring, 
the charge derealization as well as the pathways for charge 
flow. 

Summing up, after checking whether the initial requirements for 
efficient photocatalytic degradation are met, we found that BA 
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Figure 10: Schematic illustration of electro-optical processes that initiate the photocatalytic degradation of transparent pollutants on titania: 

(left) density of states, revealing the mixed character of the key MOs, (center) diagram showing the energy level alignment, the photoexcitation and 

the charge transfer, and (right) the molecular orbitals, showing the anchoring geometry and the charge localization. 



is disqualified by the poor absorption in the visible range. Of 
the other three, Ph has strongest optical absorption but it is 
outperformed by SA likely due to a more efficient electron 
transfer. We have to state that our approach cannot quantify the 
four different criteria discussed. Although the comparison with 
the experimental data is only qualitative, the present approach 
can explain some of the main features of the degradation curves 
observed. 

Conclusion 

We reported results of DFT and TD-DFT calculations 
performed on several transparent aromatic pollutants as well as 
complex systems consisting of the benzene derivatives adsorbed 
on a T1O2 nanocluster. Our goal was to answer questions such 
as: Why can colorless pollutants degrade under visible light? 
Why do some pollutants degrade faster than others? To answer 
such questions we determined the electronic structure and the 
optical spectra of the pollutant itself and found where the depro- 
tonation is more likely to take place. We optimized the geom- 
etry of pollutant-catalyst systems and shed some new light on 
the binding configurations of the benzene derivatives onto 
titania. We were able to dispel some misconceptions regarding 
the monodentate binding of the carboxy group, six-ring chelate 
binding of the salicylic acid. We demonstrated that the optimal 
binding configuration is bidentate bridging for the carboxy 
group and tridentate to adjacent titanium ions for SA. 



We explained why transparent pollutants adsorbed onto a cata- 
lyst that absorbs only in the UV can degrade under visible light 
based on the mixed character of the key MOs involved in 
optical transitions. These orbitals are delocalized on both the 
pollutant and the catalyst such that some low energy transitions 
that were forbidden become allowed. We attempted to explain 
the experimental facts which state that the efficiency of degra- 
dation under visible light irradiation decreases in the sequence 
SA > Ph > pHBA > BA. Thus, analyzing the pollutants against 
some requirements for efficient photocatalytic degradation, we 
found that BA has no absorption in the visible range even 
bound to the catalyst. Ph has best light absorption but weakest 
anchoring, pHBA has very weak electron density on the anchor, 
which hinders the charge transfer, and that SA has strongest 
binding and offers the best pathways for charge flow. 
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